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We show that accreting neutron stars in binary systems or in Landau-Thorne- 
Zytkow objects are good candidates for continous gravitational wave emission. 
Their gravitational radiation is strong enough to be detected by the next generation 
of detectors having a typical noise of 10 — 23 Hz -1 / 2 . 



1 Introduction 

A crude estimate of the gravitational luminosity of an object of mass M, mean 
radius R and internal velocities of order V can be derived from the quadrupole 
formula: 




where R s := 2GM/c 2 is the Schwarzschild radius associated with the mass M 
and s is some asymmetry factor: s = for a spherically symmetric object and 
s ~ 1 for an object whose shape is far from that of a sphere. According to 
formula (Q), the astrophysical objects for which s ~ 1, R ~ Rs and V ~ c 
may radiate a fantastic power in the form of gravitational waves: L ~ c 5 /G = 
3.6 x 10 52 W, which amounts to 10 26 times the luminosity of the Sun in the 
electromagnetic domain! 

A neutron star (hereafter NS) has a radius quite close to its Schwarzschild 
radius: R ~ 1.5 — 3 R s and its rotation velocity may reach V ~ c/2 at the 
equator, so that they are a priori valuable candidates for strong gravitational 
emission. The crucial parameter to be investigated is the asymmetry factor 
s. It is well known that a uniformly rotating body, perfectly symmetric with 
respect to its rotation axis does not emit any gravitational wave (s = 0). 
Thus in order to radiate gravitationally a NS must deviate from axisymmetry. 
Moreover, CW emission is possible only if the NS accreets angular momentum 
from an angular momentum reservoir. 

a to appear in the Proceedings of the International Conference on Gravitational Waves: 
Sources and Detectors, Cascina (Pisa), Italy — March 19-23, 1996, Eds. I. Ciufolini, F. Fide- 
caro (World Scientific, in press). 



Low Mass X-ray Binary systems (LMXB) and High Mass X-rays Binary 
systems (HMXB) are a good examples of a NS coupled with an angular momen- 
tum reservoir. These systems are formed by a NS and an ordinary companion. 
If the two stars are close enough, the NS accretes matter (and angular momen- 
tum) from the companion and consequently can radiate CW if its axisymmetry 
is broken. 

The fate of such a system depends on the mass of the companion. For the 
HMXB for which the companion is a massive star (> 8M ) and consequentely 
the life time is quite short (« 10 6 yr) the companion evolves until when the 
nuclear fuel is exhausted and becomes a supernova. If the binary system is not 
disrupted by the explosion, the outcome is a binary pulsar. PSR B1913+16 
is a good illustration of this scenario. If the companion is kicked away by the 
explosion, then the outcome is an isolated pulsar. 

If the mass of the companion is lower than 1 Mq (LMXB), the NS is 
spun up by the accretion of matter and angular momentum and, provided 
the axisymmetry is broken, the NS radiates gravitational waves steadily if the 
accreted angular momentum is evacuated via gravitational radiation. The light 
companion is evaporated by the electromagnetic emission of the NS and the 
final outcome is an isolated millisecond pulsar. 

In the intermediate case (mass of the companion between IMq and 8Af Q ) 
the final state is a binary system formed by a white dwarf and a millisecond 
pulsar. The important point is that by measuring the period modulation of 
the pulsar in a binary system it turns out to be possible to measure the mass 
of the NS. The mass of the NS is a fundamental parameter as will be explained 
later. 

Landau-Thorne-Zytkow objects (LTZO) constitute another example of a 
NS coupled_with an angular momentum reservoir. These objects, introduced 
by Landau □ to explain the stellar source of energy, have been discussed in 
details by Thorne and ZytkowH. They look as ordinary red supergiant stars, 
the main difference being the core which is a NS instead of being white-dwarf 
like. The origin of these objects (if they exist) is beleived to be a HMXB 
during the phase in which the NS is orbiting into the envelop of the companion. 
Another possible origine are aborted supernovae, i.e. supernovae for which the 
explosion is not strong enough to eject the envelop. The NS forming the core 
of these objects accrets matter from the envelop at the maximum rate, i.e. at 
the Eddington limit : Mq yr" 1 . If some amount of angular momentum is 
stored in the envelop, the NS is spun up by this accretion. The life of a LTZO 
is about 10 s yr, until the mass of the NS reaches the maximum value Af max 
and the NS collpases into a black hole. It must be noticed that the mass of 
NS, during its life, varies between the values of the mass at which the NS is 
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born to the value of M max that depends on the equation of state (EOS). 
2 Symmetry Breaking Mechanisms 

As already said, gravitational waves are radiated by a rotating NS only if its 
axisymmetry is broken. Two distinct classes of symmetry breaking mecha- 
nisms exist: The axisymmetry can be broken spontaneously (via some kind 
of instability of the NS) or the axisymmetry can be broken via some external 
mechanism. Both cases are pertinent to what follows and therefore will be 
discussed in some detail. 

A rotating NS can break spontaneously its axial symmetry if the ratio 
of the rotational kinetic energy T to the absolute value of the gravitational 
potential energy, \W\, exceeds some critical value. When the critical threshold 
Ty|W| is reached, two kinds of instabilities may drive the star into the non- 
axisymmetric state: 

1. the jChandrasekhar- Friedman- Schutz instability (hereafter CFS instabil- 
ity) a, 0,13 driven by the gravitational radiation reaction. 

2. the viscosity driven instability i. 

Let us recall some classical results from the theory of rotating Newtonian 
homogeneous bodies. It is well known that a self-gravitating incompressible 
fluid rotating rigidly at some moderate velocity takes the shape of an ax- 
isymmetric ellipsoid: the so-called Maclaurin spheroid. At the critical point 
T/|W] = 0.1375 in the Maclaurin sequence, two families of triaxial ellipsoids 
branch off: the Jacobi ellipsoids and the Dedekind ellipsoids. The former are 
triaxial ellipsoids rotating rigidly about their smallest axis with respect to an 
inertial frame, whereas the latter have a fixed triaxial figure in an ineptial 
ame, with some internal fluid circulation at constant vorticity (see ref. 13 or 
for a review of these classical results) . The Maclaurin spheroids are dynam- 
ically unstable for T/|W| > 0.2738. Thus the Jacobi/Dedekind bifurcation 
point = 0.1375 is dynamically stable. However, in presence of some 

dissipative mechanism such as viscosity or gravitational radiation (CFS in- 
stability) that breaks the circulation or angular momentum conservation, the 
bifurcation point becomes secularly unstable against the I = 2, m = 2 "bar" 
mode. Note also that a non-dissipative mechanism such as a magnetic field 
with_a component parallel to the rotation axis breaks the circulation conserva- 
tions and may generate a spontaneous symmetry breaking. If one takes into 
account only the viscosity, the growth of the bar mode leads to the deformation 
of the Maclaurin spheroid along a sequence of figures close to some Rlemann 
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S ellipsoid^] and whose final state is a Jacobi ellipsoid Eil. On the opposite, if 
the gravitational radiation reaction is taken into account but not the viscosity, 
the Maclaurin spheroid evolves close to another Riemann S sequence towards 
a Dedekind ellipsoid E3. 

The CFS instability is due to the coupling between the degrees of freedom 
of the star and gravitational waves: the star can loose angular momentum (and 
kinetical energy) via gravitational radiation. The formation of waves on the 
sea when the wind blows is due to an analougous mechanism: in the frame 
of reference of the wind, the water looses momentum because of its coupling 
with the atmosphere. Two conditions must be fulfilled to allow for the growth 
of the CFS instability: (i) the phase velocity of the pertubation must be less 
then the rotation velocity at the equator of the star, (ii) the viscosity must be 
less than a threshold value fx cr it . 

The first condition is always met: it turns out that the phase velocity of 
the gravity waves (the so-called f modes) is oc Z^ 1 / 2 , where I is the "quantum" 
number of the wave in the harmonic functions expansion (all that in a complete 
analogy with the sea waves). On the contrary, the second condition is hardly 
fulfilled: the dumping effect due to the viscosity grows as I 2 . Therefore, taking 
into the account the viscosity ^£ nuclear matter, only the mode I = 2 can 
survive. Recent computations 113 show that this kind of instability can exist 
only during a short period in the life of the star: in fact if the NS is too hot 
(resp. too cool), the bulk viscosity (resp. the shear viscosity) inhibits the 
instability. Actually the interior of a NS is more complicated: it is superfluid 
and type 2 superconductor. Superfluid vortices are coupled with magnetic 
fluxoides via their own magnetic field. Vortices and fluxoides are strongly 
pinned in the solid crust of the star. All that results in an effective viscosity 
higher than the one computed in the absence of magnetic field. Moreover, any 
mechanism that tends to rigidify the rotation of star (for example the magnetic 
field) acts against the CFS instability. From the above it turns out that the 
CFS instability seems to be very unlikely. 

The viscosity driven instability seems to be more promising: in fact, its 
rising time decreases when the viscosity increases. The physical mechanism of 
this instability is very simple: consider the rotational kinetical energy of the NS 
at fixed angular momentum L: T = L 2 / 1 where / is the moment of inertia with 
respect to the rotation axis. The kinetical energy T decreases if / increases. 

b The Riemann S family is formed by homogeneous bodies whose fluid motion can be 
decomposed into a rigid rotation about a principal axis and a uniform circulation whose 
vorticity is parallel to the rotation vector. Maclaurin, Jacobi and DcdckincUcllipsoids are 
all pppcial cases Q f Riemann S ellipsoids (for more details, cf. Chap. 7 of ref. Q or Sect. 5 of 
rcf.B). 
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It turns out that for a large enough L, the total energy of the star (sum of the 
kinetical and gravitational energy) decreases when / increases. The natural 
way to increase I is to let the configuration to be tri-axial. It is worth to note 
that the final stellar configuration is again an equilibrium configuration (in a 
rotating frame). The transition between Maclaurin and Jacobi configurations 
is a real Landau phase transition of the second order as was showed by Bertin 
and RadicatiEJ. The reader can find more details in our lecture on the subject 
at Les Houches Schoolli3. 

The main problem is that this instability can work, as already said, only if 
the NS rotates fast enough. The maximun angular velocity of a rigidly rotating 
star is achieved when the velocity at the equator is equal to the Keplerian 
velocity. It turns out that if the EOS of the fluid forming the star is too soft, the 
Keplerian velocity is less than the critical velocity for which the axisymmetry 
breaks. For a polytropic EOS (Pis. a^-japd in the Newtonian theory, 7 must 
be greater than 7 crit = 2.238 (ref.E3,E3,t3). 

3 Results for realistic equations of state 

Recently, we have generalized the-.above results to the existing "realistic" EOS 
in a General Relativistic frame E3. Table [j] shows the results: among the 12 
EOS taken under consideration, five are stiff enough to allow for the transition 
toward a 3-D configuration. In table [l[ the EOS are labeled by the following 
abbreviations: PandN refers to the pure neutron EOS of Pandharipande Ej, 
BJI to model IH of Bethe & JohnsonE], FP to the EOS of Friedman & Pand- 
haripandeB, HKP to the.no = 0.17 fm~ 3 model of Haensel et al. c3, DiazII to 
model II of Diaz AlonsoE3, Glendl. Glend2 and Glend3 to respectively the case 
1, 2, and 3 of Glendenning EOSeJ, WFF1, WFF2 and WFF3 to respectively 
the AV14 + UVII, UV14 + UVII and UV14 + TNI models of Wiringa et al. E2, 
and WGW to the A£° onn + HV model of Weber et al. El 

From the above results it appears that only NSs whose mass is larger 
than 1.64 M© meet the conditions of spontaneous symmetry breaking via the 
viscosity-driven instability. The above minimum mass is quite below the max- 
imum mass of a fast rotating NS for a stiff EOS (3.2 M© E3). Note that the 
critical period at which the instability happens (P = 1.04 ms) is not far from 
the lowest observed one (1.56 ms). The question that naturally arises is: do 
these heavy NSs exist in nature ? Only observations can give the answer; in 
fact, the numerical modelling of a supernova core and its collapse L3 cannot 
yet provide us with a reliable answer. The maximum critical rotation period 
(1.2 ms) at which the the instability appears is compatible with the rotation 
period of the fastest known pulsar (1.56 ms); moreover the age of these pulsar 
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EOS -^max -^max -Pbreak He, break -^fbreak 

[M ] [M ] [ms] [ms] [M ] 



HKP 


2 


.827 


3 


.432 





.737 


1 


.193 





.168 


1.886 


WFF2 


2 


.187 


2 


.586 





.505 





.764 





.292 


1.925 


WFF1 


2 


.123 


2. 


.528 





.476 





.728 





.270 


1.742 


WGW 


1 


.967 


2. 


.358 





.676 


1 


.042 





.170 


1.645 


Glend3 


1 


.964 


2 


,308 





.710 








stable 




FP 


1 


.960 


2. 


.314 





.508 





.630 





.412 


2.028 


DiazII 


1 


.928 


2 


.256 





.673 








stable 




BJI 


1 


.850 


2. 


.146 





.589 








stable 




WFF3 


1 


.836 


2 


,172 





.550 





.712 





.327 


1.919 


Glendl 


1 


.803 


2. 


,125 





.726 








stable 




Glend2 


1 


.777 


2. 


,087 





.758 








stable 




PandN 


1 


.657 


1. 


,928 





.489 








stable 





Table 1: Neutron star properties according to various EOS: M^jj is the maximum mass 
for static configurations, M^ x is the maximum mass for rotating stationary configurations, 
Pk is the corresponding Kcplerian period, Pbreak is the rotation period below which the 
symmetry breaking occurs, -ff Cj break is the central log-enthalpy at the bifurcation point and 
M break is the corresponding gravitational mass. The EOS are ordered by decreasing values 

of M stat 
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spans between 10 7 and 10 9 yr. 

The real problem is the minimum mass, 1.64M©, for the triaxial instability 
to develop. This is_aot in very good agreement with the measured masses (all 
in binary systems)Ej, except for PSR J1012+5307 which appears to be a heavy 
NS: 1.5 M s < M < 3.2M Q E3. Four NS masses (all in binary radio pulsars) 
are known with a precision better than 10% and they turn out to be around 
1.4 M Q El Among the X -ray binary NSs, two of them seem to have a higher 
mass: 4U 1700-37 and Vela X-l (1.8 ±0.5 M and 1.8 ±0.3 M Q respectively). 
These objects show that NSs in binary systems may have a mass larger than 
1.64 M®. 

A natural question that may arise is: why do X-ray binary NSs, which 
are believed to be the progenitors of binary radio pulsars, have a mass larger 
than the latter ones ? We have not yet any reliable answer to this question. 
A first (pessimistic) answer is that the measurements of X-ray NS masses are 
bad (compare the error bars of the masses aLthe binary radio pulsars with the 
ones of the X-ray binaries in Fig. 3 of ref. cB), and consequently not reliable. 
Actually it should be noticed that the error bars of the X-ray pulsars do not 
ve the same statistical meaning as the error bars of the binary radio pulsars 
they give only the extremum limits of NS masses in the X-ray binary. 
Consequently 1.4 Mq is not incompatible with these masses. 

A related question arises naturally: why are the observed masses of mil- 
lisecond radio pulsars almost identical ? Following the standard model, a 
millisecond radio pulsar is a recycled NS, spun up by the accretion of mass 
and angular momentum from a companion. The observed mass and angular 
velocity are those of the end of the accretion process. Consequently the ac- 
creted mass depends on the history of the system and on the nature of the 
companion. By supposing "per absurdo" that all NSs are born with the same 
mass, it is difficult to understand why the accreted mass is the same for all NSs. 
A possible answer is that this could result from some observational selection 
effect. For example, suppose that accreted matter quenches the magnetic field, 
it is then easy to imagine that the final external magnetic field depends on the 
mass of the accreted plasma. If the accreted mass is large enough, the mag- 
netic field can be lower than the critical value for which the pulsar mechanism 
works. On the contrary, if the accreted mass is quite small, the magnetic field 
is large and the life time of the radio pulsar phase is shorter and consequently 
more difficult to observe. 
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4 Detectability 



If the nuclear matter EOS is stiff enough and accreting NS in binary systems 
have a mass large enough for the symmetry breaking to take place, accreting 
NS are efficient gravitational wave emitters. It is very easy to compute the 
amplitude of the emitted gravitational waves. By equating the rate of the 
accretecLangular momentum to the rate of radiated angular momentum one 
obtains E3 

h = 1.8 x 10- (±™L) ( — \ V2 (2) 

V v J V10~ 8 erg cm- 2 s- 1 / v ' 

where h is the strain of the gravitational wave, v the rotation frequency of 
the source, Fx the X-ray flux received on Earth. Note that the distance of 
the source does not appear in the above formula. From (^), the signal-to- 
noise ratio S/N can be easily computed in terms of the observation time T 
and the sensitivity of the detector B. For the brighest X-ray source, Sco X-l 
(Fx = 2 x 1CP 7 ergem -2 s _1 ), we obtain 



Sco X-l 




(3) 



From the above formula, we see that one month of observation is sufficient to 
obtain S/N — 1 with a detector of the 10 -23 Hz -1 ' 2 class. This is however 
a misleading result: in fact, because the frequency of the CW emission is not 
known, a signal-to-noise of about 7 is required in order to have a detection 
with a confidence level equivalent to the ordinary 3c criterium. This means 
that 2.5 years of observation time with one detector are necessary to detect the 
gravitational radiation emitted by rotating NS. With 3 detectors (e.g. 2 LIGO 
+ VIRGO) the situation appears more favorable: 10 months of observation 
time instead of 30. Moreover a less naive strategy can be used to couple the 3 
detectors; we do not discuss this possibility here. 

LTZO objects are also good candidates. The radiation mechanism is anal- 
ogous to that of the accreting binary sources: the NS forming the core is spun 
up by the accreted angular momentum from the envelop. The main avantage 
of these sources is that the mass range of the inner NS spans from the initial 
mass of the NS (<~ 1.4 M Q ) up to the critical mass (> 2 M Q ). The drawback 
is that we do not if these objects exist. 

Finally, note that a deformation (ellipticity) as small as s=s 10~ 8 is sufficent 
to radiate the accreted angular momentum at the Edington mass accretion 
rate (10 -8 M Q yr _1 ). A question naturally arises: do there exist any other 
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mechanism able to deform the NS by a such a amount ? No alined magnetic 
field can do the job. The accreted matter is funelled by the magnetic field 
onto the crust of the NS, and spreads out on the surface, but magnetic field 
acts as a magnetic brake for this process. The efficiency of this magnetic 
braking depends on the conductivity of the plasma and on the strength of 
the magnetic field. A rough estimation of the typical spreading time r of 
the accreted plasma on the surface of the NS gives r ^> 1 yr. This means 
that 3-D asymetries can be larger than 10~ 8 for the accreting rate of 10~ 9 — 
1CP 8 M Q yr _1 . The above encouraging result is correct only if the effecive 
conductivity of the plasma is equal to the microscopic one. Indeed plasma 
instabilities can reduce the effective conductivity by orders of magnitude. The 
most dangerous of them is the instability generating the reconnection of the 
magnetic field lines. Fortunately, no X or O point exists in the magnetic field 
configuration. Therefore this kind of instability seems to be excluded. More 
investigation is needed to clarify this important question (the Authors thank 
Prof. E. Spiegel and Dr. A. Mangeney for illuminating discussions on this 
point). 

5 CONCLUSION 

Accreting NS in binary systems or in LTZOs can be good gravitational CW 
emitters. Their positions on the sky are known, therefore data can be easily 
reduced to the solar system barycentric frame and the Doppler shift induced 
by the motion of the Earth can be properly taken into account. The amplitude 
of the predicted gravitational waves is large enough to be detected with the 
1CP 23 Hz -1 / 2 class of detectors. Positive detection will give to us important 
informations on the equation of state of the nuclear matter in NS. The proof 
of existence of the LTZOs will be a major discovery leading to important 
informations on the stellar evolution during the common envelop phase. 
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